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Quantum information protocols utilizing atomic ensembles require preparation of a coherent spin 
state (CSS) of the ensemble as an important starting point. We investigate the magneto-optical 
resonance method for characterizing a spin state of cesium atoms in a paraffin coated vapor cell. 
Atoms in a constant magnetic field are subject to an off-resonant laser beam and an RF magnetic 
field. The spectrum of the Zeeman sub-levels, in particular the weak quadratic Zeeman effect, 
enables us to measure the spin orientation, the number of atoms, and the transverse spin coherence 
time. Notably the use of 894nm pumping light on the Dl-line, ensuring the state F = 4, rar = 4 
to be a dark state, helps us to achieve spin orientation of better than 98%. Hence we can establish 
a CSS with high accuracy which is critical for the analysis of the entangled states of atoms. 

PACS numbers: 32.30.Dx, 32.80.Bx, 76.70.Hb 



I. INTRODUCTION 

The spin state of ensembles of alkaline atoms have 
been studied for a long time in many different contexts, 
e.g. sensitive magnetic field measurements pj, frequency 
standards 0, and recently within the field of quantum 
information. Coherent transfer of states of the electro- 
magnetic field to atomic spin states with the aid of elec- 
tromagnetically induced transparency has been demon- 
strated 0,0]. A sample of cesium atoms has been spin 
squeezed 0, two samples of cesium atoms have been 
entangled ^3], and in 0, Q the sensitivity of atomic 
spin states to quantum fluctuations of the electromag- 
netic field was demonstrated. Some of these experiments 
require the preparation of a coherent spin state, CSS, 
(e.g. all atoms pumped into the magnetic substate F = 4, 
rriF = 4 in the cesium ground state). 

This paper provides a detailed report on the creation 
and characterization of the coherent spin state used in the 
recent experiments on entanglement and quantum mem- 
ory 0,0 13| The importance of the CSS of a macroscopic 
ensemble is primarily in that the transverse spin compo- 
nents of this state are in the minimum uncertainty state 
and thus provide a necessary starting point for observing 
quantum effects. For example, in order to apply the nec- 
essary and sufficient conditions of the entanglement and 
spin squeezing 0, 0] one has to know the variance of 
the CSS. Unlike in case of the shot noise of light, atomic 
coherent state noise cannot be easily established via inde- 
pendent measurements. It therefore has to be determined 
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experimentally in each case. 

In this paper we report measurements of the orienta- 
tion, the coherence time, and the number of atoms in 
a spin ensemble in one selected hyperfine ground state 
(e.g. the F = 4 states in cesium). This is done by 
inducing weak transitions among the Zeeman sublevels 
by radio-frequency magnetic fields. Other similar prob- 
ing methods have been demonstrated in where mi- 
crowave magnetic fields on the hyperfine transition are 
used to characterize the efficiency of optical pumping 
into e.g. the F = 3,tof = state. Also, fluorescence 
from atoms excited by the pumping process can give in- 
formation about the excited spin state which again hints 
on the ground spin state In the literature, the spin 
state of atoms is often modelled by numerical solutions 
to rate equations. We make a different approach by tai- 
loring simple ad hoc models to describe the main features 
of the spin state. 

The paper is arranged as follows. In section [H] we de- 
scribe the ground spin state of our atoms theoretically. 
This includes an introduction to the notation and the 
particular physical system used in the experiment III Al 
deriving equations of motion III Bl and modelling the dis- 
tribution and coherence time of atoms among Zeeman 
suhlevels lll CI In the experimental section liTfl we describe 
the actual setup IhTaI discuss the conditions for resolv- 
ing the quadratic Zeeman effect IIII Bl and demonstrate 
that our models actually describe the experiments IIII (11 
Much of the emphasis is put on the quadratic Zeeman 
effect and the ability to resolve this spectroscopically. 
However, in some experimental conditions this is not the 
case and we discuss how we can employ our techniques 
in a non- resolved regime IIII PI and also in a regime with 
pulsed lasers IW1 In appendix El we review the quadratic 
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as the quadratic Zeeman effect, for completeness we have 
summarized the most important results in appendix El 

Applying RF-magnetic fields we can induce transitions 
between magnetic sublevels. Our spectral resolution is 
sufficiently high in order to observe the small energy 
shifts caused by the quadratic Zeeman effect, see Fig. 
We see eight distinct lines corresponding to transitions 
between two adjacent levels among the 2F + 1 possible 
levels in the F = 4 state. The connection between the 
spin state and the spectrum in Fig. will be derived in 
the following, we will be able to fully characterize the 
spin state in many situations important for other exper- 
iments. 

We characterize the state of atoms by the density op- 
erator pij given by 



FIG. 1: An experimental spectrum (dots) of magnetic tran- 
sitions among the nine sublevels of the F = 4 ground state 
in cesium. An external RF-magnetic field (the frequency of 
which is shown on the abscissa) modulates the spin state. The 
recorded response, called the magneto-optical resonance sig- 
nal (MORS), is shown on the ordinate axis. The solid line 
is a fit to a model to be developed in section III CI where the 
parameters and the definition of MORS will also be carefully 
explained. The many peaks tells us that atoms are distributed 
among all nine levels resulting in a low orientation p — 0.346. 
The line width 9.4Hz is a direct measure of the decay rate 
of spin coherence. According to Eq. l|A5|l the corresponding 
F = 3 signal is approximately 1kHz away and does not inter- 
fere here. The inset shows the level structure of the ground 
states of cesium following Eq. <A2fl . Our experiment is carried 
out around 10 _4 T which is far into the linear regime. How- 
ever, with a sufficiently good resolution the quadratic effect 
is visible. 



Zeeman effect. 



II. THEORETICAL DESCRIPTION 
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where i,j = —F, —F + 1, . . . , F parametrizes the mag- 
netic sublevels \i) or \j) in Dirac notation, and the sum 
is done over all individual atoms. With the z-axis as 
quantization axis we may express the total macroscopic 
angular momentum of the atoms in the hyperfine state 
F as 
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In this section, we introduce the physics and notation 
which will enable us to understand how magneto-optical 
resonance can be used to characterize the spin state of 
an atomic gas sample. 



A. Spin State Evolution in the Process of the 
Double Magnetic Resonance 

Let us consider an ensemble of atoms in a Zeeman mul- 
tiplet in the ground state of an alkali, in our case this is 
the F — 4 state in cesium. Applying an external mag- 
netic field, the energy of each magnetic sublevel changes 
as shown in the inset of Fig.^ For low fields the splitting 
between two levels grows linearly with field strength with 
small quadratic corrections given by ^qz = 2^/ i/ hfs I 
where is the Zeeman splitting, or the Larmor fre- 
quency, and i^hfs is the hyperfine splitting. This is known 



where C (F, m) — y/ F(F + 1) — m(m + 1) and F± are 
raising/lowering operators for the spin along z. Note, if 
we are interested in calculating the spin components of 
F, it suffices to calculate the diagonal terms and the first 
off-diagonal terms in the density matrix. This enables us 
to describe effects like the orientation of the spin state, 
see Sec. Ill CI In general, to describe aspects of alignment 
we also need second off-diagonal terms. We restrict our 
description to diagonal and first off-diagonal terms in the 
following. 

The Hamiltonian of spins subject to a magnetic field 
B can be written 

H = g F fi B F -B + 0{B 2 ) (3) 

where the second order correction is described in ap- 
pendix El We will partly include this second order cor- 
rection, since it will prove to be important for the energy 
levels. We will place a constant bias magnetic field with 
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strength -Bbias along the z-axis and apply an RF-magnetic 
field \B RF \ cos(ui+0) along the x-axis. The Hamiltonian 
may now be written (for a single atom) 



r in in 
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where B RF — \B RF \e~ l<t> is the complex amplitude. The 
first term is primarily F Z B Z , but we take the second order 
corrections into account by explicitly stating the energy 
levels hui m of the m'th sublevel. The second term is F X B X 
originating from the x-polarized RF-magnetic field which 
induces transitions between the magnetic sublevels. The 
rotating wave approximation has been made here. We 
may wish to write the Hamiltonian entirely in terms of 
the density operators p^: 



where we have just inserted the Hamiltonian © into © 
and added the decay term, — T/2 ■ p\ 2 . We will restrict 
ourselves to a description of spins in the case where 
F x ,F y -C F z , i.e. the angle 9 that the spins deviate 
from being oriented along the z-axis is much less than 
unity. From equations Il2al2cl) this can roughly be writ- 
ten as p m +i,m ~ ■ p m ,m, and following the same lines 
Pm+2,m ~ 9 2 ■ Pm.m- It is then justified to neglect the co- 
herences P02 and P13 in the above equation. For brevity 
we will define uj 21 = o-i 2 —o-ii, which is the frequency corre- 
sponding to the transition from m F = 2 to m F = 1. This 
frequency is the Larmor frequency and it is fast compared 
to the inverse time scale for dynamical evolution of the 
spin state. Since the RF-magnetic field frequency co will 
be in the vicinity of oj 2 i it will be convenient to define 
the slowly varying operators 
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The equations of motion are now determined by 
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where the first term is the coherent evolution of the sys- 
tem, and the interaction with the environment will be 
put in by hand as decay terms. 



B. Solution of equations of motion 

We will now solve equations (0 and © , and to illumi- 
nate the method for solving these equations, we will pick 
out a single example and work it out thoroughly. The 
time derivative of e.g. pi 2 is 
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Using this definition, equation {JJ will turn into 

~ f>12 - = (i[w- ui2i] - r/2)/5i2 

"™'-C(F, l)B RF [p 22 — pu] 



dt 
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The constant T will describe the decay of the transverse 
spin components. We will assume the time scale T 2 for 
the transverse spin decay is much shorter than the time 
scale Ti for the longitudinal spin component (along the z- 
axis). Experimentally we typically find Xi » 200-300ms 
and T 2 < 40ms and the approximation is justified. Then 
the operator p~i 2 will follow (p 22 —pn) adiabatically which 
is expressed mathematically by setting dp\ 2 /dt = in 
the above equation. Alternatively, if constant pumping 
maintains F z the small angle condition 9 <C 1 will ensure 
constant {p 22 — pn). Then the above is simply the steady 
state condition where transients have been damped away. 
The result for pi 2 is then 
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This method applies to all density operators p m ,m+i, and 
substituting into equations l!2at and Il2bj) we get 
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These equations can be interpreted as a number (2F) of 
two-level systems that all interact with the driving RF- 
magnetic field. Two adjacent magnetic sublevels to and 
m+1 act as one two-level atom with the usual Lorentzian 
response (resonance frequency w m +i, m and line width 
r m+ i >m FWHM). Each two-level system does not re- 
spond with exactly the same weight which is reflected 
in the factor F(F + 1) — to(to + 1). All the resonances 
add up coherently to give the overall response of the spin 
state to the RF-magnetic field. Note, that F x and F y 
oscillate at the driving frequency w and not the "natu- 
ral" frequencies w m +i, m . This is a steady state behavior, 
just like the forced harmonic oscillator when transients 
have been damped. In section IIVI we will comment on 
non-steady state behavior of the spin system. 



C. Modeling the spin state 

In the previous subsection we derived how the spin F 
responds to an external RF-magnetic field Brf- Our 
main motivation is to use this knowledge to characterize 
the spin state, i.e. we wish to apply the field Brf and 
measure the response F x or F y in order to deduce param- 
eters like p m , m , r m+ i m and so on. Now, for cesium in 
the F — 4 hyperfine ground state there are nine popula- 
tions pm.m and eight line widths r„ l+ i jTO together with 
the resonance frequencies. To fit an experimentally mea- 
sured spectrum (see e.g. Fig. ^| to all these parameters 
will be very hard and in the following we will develop a 
model to significantly reduce the number of free parame- 
ters. We will just tailor a model and the justification for 
this model will be an experimental test. 

Let us consider a case where we wish to orient all 
atomic spins along the z-direction, i.e. attempt to put 
many atoms into the to = F substate. This can be done 
experimentally by illuminating the atoms with circularly 
polarized light, as will be described in the experimental 
section. It is then convenient to define the orientation p 
as an order of merit 



ife (12) 
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Note, that with this definition p = 1 if all atoms are in 
the extreme to — F sublevel, and p = for a completely 
unpolarized sample with p m ,m = l/(2F + 1) for all to. 
We try to let this be the only parameter describing the 
relationship between the nine populations p m . m . With 
the condition p m . m = 1 we have thus reduced eight 
free parameters to a single one. 

Now, we describe ensembles of atoms and given p we 
will assume that the spin state maximizes the entropy 
S = — Tr(/51np). To find the individual p m . m we use the 



method of Lagrange multipliers. We must solve 
d 



dp m ,r, 



S - a p m , m - /?y^m • p m ,7n) = 



(13) 



We now just need to adjust a and in order that 
Tr(p) = 1 and p is as desired. Doing this is more or 
less a computational problem and in principle not dif- 
ficult. For the physical understanding we just need to 
remember that we can write p m ,m oc e m where e = e - ' 3 is 
a parameter which is a function of p. This can go directly 
into equations Ijllafl and ijllbjl . We stress that the above 
considerations are ment as an intuitive guideline to un- 
derstand what is really just an experimentally justified 
model. This model is only valid in some specific cases, 
our example with circular polarized pumping light. Our 
experience tells us that the model works fine for large 
longitudinal relaxation times T\ and not too strong in- 
fluence by the probe laser. 

For the eight line widths r m+ i m in the case of cesium 
we will make a model with two free parameters. First, a 
common line width r com is assigned to all transitions in- 
dependent of to. The physical cause for this type of decay 
could be magnetic field inhomogeneities, collisions, and 
loss mechanisms common to all atoms. In addition, if we 
wish to create a well oriented sample with to approach- 
ing F we will need to illuminate the atoms with resonant 
circularly polarized light. In our case this is the 894nm 
6S1/27 F — 4 to &P\/2, F 1 — 4 line. This light causes ex- 
citations from the atomic ground sublevel to with a rate 
7 m cx | (4, to, 1, 1|4, to + 1) | 2 = (4-m)(5 + m)/40, where 
the second term is the square of Clebsch-Gordan coeffi- 
cients. For a magnetic transition between ground sub- 
levels to and to + 1 the resonant pumping light will con- 
tribute to the line broadening proportional to 7m+7m+i- 



Thus we define the width T p 
pumping process such that 



caused by the optical 
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pump 



(14) 



where the normalization is such that for the to = 3 — > 
to = 4 transition we have T4 3 — T com + r pump . 

Finally, we must have the resonant frequencies as pa- 
rameters in our model. We will write this as a central 
frequency w cen ter and a splitting U! sp in such that 



^m + l,m — ^center T" ^split 



(15) 



Theoretically we should have w sp iit = 2tt ■ vqz (see equa- 
tion l|A7J) l but we keep it as a free parameter since in 
practice this splitting may differ slightly from the theo- 
retical value because of Stark splittings. 

To sum up, a possible description of the ground spin 
state involves the total spin F z and the orientation p 
together with the line widths r com and r pump , and the 
frequencies w cen ter and w sp i; t . An equivalent but compu- 
tationally easier way to represent F z and p is to use the 
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FIG. 2: (a) The lasers used in the experiment together with the level scheme of cesium. The F = 4 ground state is with 
its nine sublevels of main interest in this experiment, (b) The experimental setup. The atomic ensemble is illuminated with 
optical pumping and repumping lasers with <r+-polarized light parallel to a static bias magnetic field. This orients the spins 
along the z-axis and gives rise to a 325kHz Larmor precession. An RF-magnetic field perpendicular to the spins induce small 
spin components in the plane orthogonal to z which are read out by polarization rotation of a linearly polarized probe laser 
beam (giving an AC-signal). The spin along z can also be read out on a probe laser beam (giving a DC-Faraday signal). 



number of atoms N4 = Np^ of atoms in m = 4 as one 
parameter and the parameter e such that the population 
N m can be expressed as N m = Np m , m = N^e^™^ . 



III. EXPERIMENT 

In this experimental section we will first describe the 
experimental setup in order to understand what is actu- 
ally measured and how the measurements relate to the 
equations of section [H] Having this in place we comment 
on the requirements that actually allow us to resolve the 
individual sublevels as seen in Fig. 



A. Experimental setup 

The experimental setup and the laser level scheme is 
shown in Fig. |2j Referring to part (a) of the figure, we 
have a Ti:Sapphire laser probe detuned by A w 1GHz 
from the &S\/2,F = 4 — » 6P 3 /2,F — 5 transition at 
852nm. This laser is sensitive to the state of atoms in 
the F = 4 ground state and will be used to probe these 
atoms. To control the state of atoms in the F — 4 ground 
state we have two home built grating stabilized diode 
lasers. The repump laser is tuned to the 6S1/2, F = 3 — > 
6P3/2,F — 4 transition at 852nm. This mainly serves 
to remove atoms from the otherwise dark F = 3 ground 
state. The optical pump laser tuned to the &Si/ 2 ,F — 
4 — > QPi/ 2 ,F — 4 transition at 894nm plays the most 
important role in distributing atoms among the sublevels 
of the F — 4 ground state. Note, that atoms in the 
extreme F = 4, mp = 4 state will not absorb light from 
any of the two lasers. 

Now turning to Fig.JSb, we place a paraffin coated va- 
por cell containing cesium inside a shield of ^-metal to 
isolate the atoms from laboratory and Earth magnetic 



fields. A static magnetic field of strength approximately 
0.9 Gauss is then applied in the z-direction and the re- 
pump and optical pump lasers are illuminating atoms 
along that direction. The polarization of these lasers can 
be adjusted, and by choosing circular polarization we can 
orient the atomic spins along the z-direction. 

The probe laser is split into two beams. One of them 
is directed along the x-direction transverse to the atomic 
magnetization created by the repump and optical pump 
lasers. This probe laser is linearly polarized along the z- 
direction and will undergo polarization rotation propor- 
tional to the spin F x along the x-axis (see e.g. jsllj])- 
By suitable arrangement of a A/2-plate and a PBS, the 
difference signal of two photo detectors will be propor- 
tional to the polarization rotation (for small angles). 

The atomic spins F x along the x-axis will have zero 
mean value unless an RF-magnetic field is applied trans- 
versely to the static magnetic field. We apply such a 
field at frequency uj (we shall call this frequency the lo- 
cal oscillator frequency) , and the motion of F x will now 
be described by equation (1 1 1 all . We may write the out- 
coming AC-signal as i(t) = a- F x — a- Re{A(£)} where a 
is a constant depending on beam geometry, laser detun- 
ing and intensity 0], and A(t) reflects the curly bracket 
of equation UllaJ) . We know from this equation that 
A{t) = A{uS)e lut will posses only a single frequency com- 
ponent, namely the local oscillator frequency u driving 
the transverse spins F x and F y away from zero. The 
amplitude of this frequency component is experimentally 
measured by inserting the photo current i{t) into a lock- 
in amplifier and decomposing the signal into sine and 
cosine components: 

i{t) = a ■ Re{A(u;)e luJt } 
= a ■ (Re{A(to)} cos(wi) - lm{A(uj)} sin(wi)) (16) 

The lock-amplifier may give the sum of the squared am- 
plitudes of the sine and cosine components which in our 
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case will be exactly a 2 |A(u;)| 2 . We shall call this sig- 
nal our magneto-optical resonance signal at frequency u> 



J 



MORS(w) = const 



F-l 



i(uj 



(MORS(w) in short). Combining the above with equa- 
tion 111 la! and ignoring irrelevant constants we find 



[F(F + 1) - m(m + 1)] 



Vn+l,m 



u) 



m+1 



[Pm+Lm+1 



(17) 



The second part of the probe laser beam in Fig. 03 is 
linearly polarized and is directed along the static mag- 
netic field direction and will be subject to polarization 
rotation proportional to the atomic spin F z . The angle 
of rotation 6*dc (denoted the DC-Faraday rotation angle) 
can be measured and we have in analogy with II17J1 



#dc = const • F z = const • N 



F 



TTlpm,Tn 



(18) 
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With the magneto-optical resonance signal, the DC- 
Faraday rotation angle, and the methods of section ITTCl 
at hand we will be able to say much about the number 
of atoms N of the sample, the populations p m ,m and the 
decay rates T rn +i. m of the coherences p m ,m+i- 



B. Resolving the different Zeeman lines 

To resolve the different Zeeman lines we obviously need 

Fm+l,m 

< cjqz- There are many contributions to the 
decay of the transverse coherences, collisions between 
atoms, collisions with the cell walls, power broadening by 
laser light, dephasing by magnetic field inhomogeneities. 

Atoms are kept in a paraffin coated glass cell which 
prevents atoms from depolarizing when hitting the walls. 
This method can increase the coherence time up to 
roughly 1 second 0]. However, our record is somewhat 
lower at 50ms. 

The probe laser beams give only a marginal contribu- 
tion to the decay if they are detuned sufficiently far from 
the resonance or has a low intensity. We estimate the rate 
r p h at which an atom scatters photons by the two-level 
atom result (see e.g. [lfl| l 
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(19) 



i+{2\/- i )' 1 ls t ne saturation parameter. / is 
the beam intensity with 7 sa t = 2Tr 2 hc^/3X 3 being the 
saturation intensity. A is the optical wave length, 7 is 
the natural line width of the optical transition, and A is 
the detuning (assumed much greater than 7 in the last 
step of the equation). In our case we operate at up to 
lmW/cm 2 and the detuning is typically around 1GHz 
which contributes the order of T p h = 100Hz to the decay 
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FIG. 3: The measured line width as a function of an applied 
magnetic field gradient. We confirm the scaling derived in 
equation 12(H . The 8.7Hz minimum is set by collisions, power 
broadening or possible residual inhomogeneities and the ad- 
ditional quadratic part stems from the applied magnetic field 
inhomogeneity. 



rate according to the above estimate. We usually ob- 
serve a somewhat smaller width which probably can be 
attributed to the crude two-level atom approximation. 
Nonetheless, to be able to see the magneto-optical res- 
onance signal a considerably weaker probe can be used 
and the power broadening can easily be reduced below 
1Hz. 

The repump laser beam is detuned roughly 9GHz from 
the F = 4 ground state and thus plays a very minor role 
even at strong intensities. The optical pump laser is di- 
rectly on resonance with the F = 4 ground state sublevels 
and then only a moderate intensity can be allowed if dif- 
ferent Zeeman lines should be resolved. The quantitative 
aspects of the optical pump laser power broadening was 
discussed in equation 11411 . 

The magnetic field must be homogeneous to a high 
degree. Since atoms are moving around in the magnetic 
field, different atoms will experience different magnetic 
field strengths and hence different Larmor frequencies. 
To estimate the role of a possible gradient dB/dz we use 
the following simple model (see e.g. 16]). First, divide 
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the atomic sample into two parts, 1 and 2, along the bias 
magnetic field direction. If the sample length is L and the 
bias field has strength B , the field strength in the two 
parts will be of order Bq±OB jdz ■ L and the difference in 
Larmor frequency will be {gF^B /h)8B / 'dz ■ L according 
to l|A4l) . We follow an atom during the time T it takes 
for it to decohere. If v is a typical speed of the atomic 
motion, the number of visits ni in part 1 or n 2 in part 
2 will be of order Tv/L, since each visit has duration 
L/v. The difference n\ —n 2 has mean zero and standard 
deviation of the order std(ni — n 2 ) = yTvjL. Thus the 
uncertainty 5<p in the accumulated phase during Larmor 
precession is 

xj.~9fVb dB L 

— h — ' ~dz ' ~v ' ' ni _ n2 > 

In the first line we set 5<j) equal to unity since this is the 
situation after the time of decoherence T. We see that 
the broadening T^h by inhomogeneities scales quadrat- 
ically with the field gradient. If we take gF ~ 1/4 (see 
eq. L — 0.030m, v = \JksT /tocs = 137m/s at 

T = 300K, we get <?fA*b/^ = 350Hz/mG and expect the 
broadening to be T iah = 0.024Hz • m 2 /mG 2 • (dB/dz) 2 . 

The experimental investigation can be seen in Fig. 
and we definitely confirm the scaling law predicted above. 
The numbers match within a factor of two which puts 
some confidence to our simple model but this is probably 
also partly luck since we were very crude in the model 
w.r.t. factors of 2 and n. Comparing the experimental 
result with the splitting due to the quadratic Zeeman 
effect (G3> we find for our particular setup that in order 
to have r^h < vqz we must have 1/B ■ dB/dz ■ L < 
1.2 -10" 3 . 



C. Confirming the spin modeling 

We will now give experimental support to the theo- 
retical derivations and the modeling of the spin state as 
described in section [H] We first focus on the validity of 
the simple model introduced in subsection III CI i.e. we 
test the dependence p m ,m oc e m , the model for the line 
widths irnfji . and the frequency distribution l(T5|l . 

With the setting as in Fig. 03 we leave the optical 
pump laser off and vary the polarization of the repump 
laser. This will create different distributions of p m ,m and 
we record the MORS. In Fig. we see an example where 
all eight lines are clearly visible, the polarization of the 
repump laser was here quite far away from being purely 
circular. The dots are experimental points and the solid 
line is a fit to the model l|17fl with N, e, T com , w cen ter, and 
w sp iit as free parameters. r pump is set to zero. We see 
that the solid line matches the experimental points very 
well. Note, that e corresponding to p = 0.346 is the only 
parameter really describing the relative strength of the 



individual peaks, while the other parameters are common 
to all peaks. The center frequency w con t or « 325250Hz is 
set by the magnetic field (or to be true, this frequency was 
a quiet place in terms of laser noise for other experiments 
and thus our detectors was tuned to this frequency) . The 
splitting w sp ut = 22Hz is close to the 23Hz expected from 
Eq. (|A7j) . The small deviation is due to Stark shifts from 
the laser beams. Finally, we find r com = 9.4Hz (FWHM). 
This corresponds to a life time of the spin coherence of 
34ms. 

The next example is recorded with the repump laser 
purely circular and optical pump laser still off. The spec- 
trum can be seen in Fig. Now the spectrum is dis- 
placed much more to one side and the fit gives p = 0.823. 
This single parameter still seems to describe the shape 
with good accuracy. 

The third example we will show is seen in Fig. QJa. 
Here the situation is as before but now with a weak opti- 
cal pump present. We observe an additional broadening 
of the left most peak by r pump = 5.5Hz and we also 
note that the second peak seems much broader (should 
have an additional broadening by 15.1 Hz according to 
Eq. ITljl ). Since the fit and the experimental points fol- 
low each other very well, we get support for the modeling 
of r pump . The orientation p — 0.967 shows that we are 
very close to have all atoms in m_p =4 with only a mod- 
erate amount of optical pumping light. 

In the three examples described above we also get a fit 
for the number of atoms N on a relative scale (the con- 
stant in front of Eq. I|17J1 depends on many experimental 
parameters so we do not wish to calculate it in absolute 
units). Multiplying this N with the fitted orientation p 
gives the total spin J z (on a relative scale). Now, the 
DC-Faraday rotation signal 6*dc (see iflSjl l gives an inde- 
pendent measure of J z and we may compare the fitted J z 
with 6*dc to get another consistency check on the model. 
This is shown in Fig. where we plot the fitted J z as a 
function of One- The lowest points are taken with the re- 
pump laser only and varying repump polarization. The 
upper eight points are taken with purely circular opti- 
cal pump of increasing intensity in addition to a purely 
circular repump laser. We see a very nice agreement be- 
tween the fitted and the directly measured values giving 
strong support to both the derivations leading to Eq. Q17J1 
and the modeling of the spin state described in subsec- 
tion mo 



D. Unresolved lines 

The spectra shown so far have been more or less well 
resolved which enabled us to directly determine the ori- 
entation p. Now, how much information can we extract 
if the line widths are much broader than the quadratic 
Zeeman splitting w sp iit? 

First, assume that all atoms are subject to decoherence 
with the same rate described by T com 3> w sp ut and decay 
from pumping light is a small contribution r pump w 0. 
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FIG. 4: Two examples of experimental (dots) and fitted (solid line) traces. The left graph was obtained with pure cr-polarized 
repump laser and no optical pumping. On the right graph a small amount of optical pumping light is added giving rise to a 



non-zero F D 



Here one can, on careful inspection, observe that the two leftmost peaks have different widths. Note also, that 



the height has grown by a factor of three compared to the graph on the left. 




FIG. 5: (a) Consistency check of the models. Fits to different spectra give an estimate of J z . This can be compared directly 
to the independently measured DC-Faraday rotation signal foe which is proportional to J z . We indeed observe a straight 
line through the origin. Note, the model description p m ,m <x e m proved to be less accurate for the lowest four points, but by 
coincidence the points still fit well, (b) The interdependence of p, J z and r pump in the limit where r pump dominates both the 
common width F com and the quadratic Zeeman splitting w sp ut. Different choices of fitting parameters will lead to a satisfactory 
fit and additional information is needed to place precise bonds on p. If for example J z is measured independently with an 
accuracy of 2% the orientation can be defined within 1% in the example shown. 



Then Eq. JTJJ reduces to 



MORSH = const' 



2N 



J2 m Pm,r 



i(uj - to) + r c 



« \J Z \ 2 (21) 



We see that in this case the spectrum will be a single 
Lorentzian the size of which is only depending on J z . In 
this case the independent measure from the DC-Faraday 
signal would only contribute exactly the same informa- 
tion and we will not be able to deduce the orientation 
P- 

On the other hand, if r pump dominates r com and 
^spiit we will get a signal that depends on the internal 



atomic spin state. To examine this approximation we set 
Tcom = w S piit = and try to fit the rest of the param- 
eters to a spectrum which is a perfect Lorentzian. The 
correct fitting parameters of course have p — 1 and r pU mp 
equal to the Lorentzian width but in practical life other 
sets of parameters will also fit the spectrum to some ex- 
tent. We find that orientations in the range p = 0.9 to 
p = 1 result in an agreement one would find reasonable 
if the spectrum was an experimental trace. By fixing p 
to a value in this range, the values given from the fit 
of J z and r opt are shown in Fig. 0D. We see that if we 
can estimate one of the parameters J z or r op t indepen- 
dently we should be able to calculate the orientation p. 
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FIG. 6: An example of the magneto-optical resonance signal 
in a pulsed regime. The spin state is close to the maximal 
state with uif = 4 as assumed in the text. Note the ripples 
which have frequency spacing (15ms) _1 . The decay rate F 
takes different values during one measurement cycle depend- 
ing on which lasers are on. The timing of laser pulses is shown 
in the inset and the corresponding decay rates shown on the 
left part of the figure. 



For instance, a measurement of J z to an accuracy of 2% 
will fix the orientation p to 1%. One only needs to have 
one fix point, e.g. if one knows that we have p = 1 per- 
fectly in one case, or if one can reduce r pump to the point 
where the spectral lines become resolved and a calibra- 
tion like Fig. can be performed. Experimentally we 
have seen orientations better than p = 0.98 for atomic 
densities around 10 11 cm~ 3 . The optical pumping laser 
at the 894nm Dl-line is essential to this achievement. 
We have tried to optically pump on the D2-line with 
somewhat lower orientation as a result (a little above 
p = 0.9). A possible explanation is that the rescattered 
light on the F — 4, mp = 4 — > F = 5, mp = 5 tran- 
sition from one atom affects the state of other atoms. 
Indeed, according to [ijj, even with a dark state when 
using 894nm pumping light one would expect problems 
with densities higher than pc — (fxi?) -1 because radia- 
tion will be trapped inside the sample. Here a is the cross 
section for light absorption and R is the extent of the gas 
sample. Our atomic sample is Doppler broadened with 
the width Spd = 378MHz. With a natural line width of 
the 894nm Z?l-transition of 7 = 4.6MHz and a sample 
extent of R = 3cm we estimate the critical density pc to 
be roughly pc ~ [\ 2 /2n-j/5v B - R}- 1 = 2-10 n cm- 3 . We 
see that we are in the regime where radiation trapping 
may be a limiting factor, but the experiments tell us that 
the limitations are still small. 



IV. PULSED EXPERIMENTS 

All previous derivations and measurements have been 
carried out in cw settings, i.e. Eqs. UllaJ) and i|llhjl as- 
sume constant values of frequency u> and decay rates 



r?n+i,ra- This is indeed valid if lasers are running cw 
and if we scan the frequency oj slowly enough. But some 
experiments must be carried out in a setting with pulsed 
lasers, e.g. one might wish to prepare the spin state in 
the maximally oriented state F — 4, mp — 4 by illumi- 
nating atoms by a pulse of resonant, circularly polarized 
laser light. Such state preparation has been used, e.g. for 
creation of entanglement between two samples of cesium 
atoms 0. For the magneto-optical resonance method to 
be useful in such experiments it must be utilized in the 
correct experimental conditions which now means time 
varying decay rates T m +i tm . In this section we outline 
the extensions into the pulsed regimes and discuss the ap- 
plicability of the magneto-optical resonance method for 
characterization of spin states under these circumstances. 

We assume for simplicity that atoms are pumped to 
mp = 4 to an extent that we only need to consider tran- 
sitions between mp — 3 and mp = 4. The extension to 
all levels should be straightforward (but cumbersome). 
For these two levels we may write Eq. (jHJ as 



dm 

dt 



(iA - r/2) j 5 34 + ix[p4A - P33] 



(22) 



where A = ui — W43 and \ — gppsBpipC(F,3)/ Ah. We 
assume as in section III Bl that the populations p44 and 
P33 can be treated as constants corresponding to small 
angle deviations from the z-axis. Then the solution of 
the above equation is straightforward 



P34(i)=/534(0)e^- r / 2 ) t 

-JA^^-^-^-^) 



(23) 



This solution starts out with ^34(0) at t = and 
makes a damped oscillation toward the steady state value 
—ix[p44: — P33]/(*A — r/2). Note, this steady state value 
is exactly the result in ifTojl . and it is reached in a time 
w r _1 . With the solution of P34 we can continue to find 
the actual spin, e.g. F x given by Il2at and predict the 
results of a measurement. 

Experimentally, we set up pumping lasers and a probe 
laser measuring the transverse spin state as in Fig. Eh- 
The lasers are turned on and off with acousto- and 
electro-optical modulators. The decay rate in the absence 
of lasers is denoted Tdark which is typically small. When 
the probe laser is on an additional broadening of r pro b e 
is present leading to a total decay rate of r pro b c + Tdark- 
During the optical pumping pulse the atoms are typ- 
ically subject to a high decay rate given in total by 
r pU mp + Tdark- The probe laser is typically turned on 
shortly after the optical pumping has been turned off 
and is maintained for a time shorter than the decay time 
(r pr obe + Tdark) -1 - We are thus in the transient regime 
of Eq. 123t and given the frequency u> of the driving RF- 
magnetic field we cannot obtain a simple estimate of the 
amplitude of the response at that frequency as in itTTjl . 
Instead we have time varying quadrature components of 
the measured AC-signal and we simply average these over 
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time in the presence of the probe laser. From the perspec- 
tive of modeling we need to evolve ^34 according to ffify 
with the relevant decay rates and integrate the result over 
the time of the probe laser pulse. 

An example of an experimental trace with correspond- 
ing numerical modeling is shown in Fig. U3 The inset 
shows the pulse sequence of lasers with a total period of 
15ms. The repump and optical pump lasers are turned 
on for 1.0ms with a high power making a power broad- 
ening r pump = 770 ± 30Hz. After a 0.3ms period with 
lasers off the probe laser is fired for a 0.5ms period giv- 
ing rise to a small power broadening r pro b ~ 2Hz in 
addition to the background value around Tdark ~ 18Hz. 
These decay rates are estimated from cw experiments car- 
ried out with the same laser powers. The experimental 
trace matches quite well the calculated spectrum in the 
sense that the general trend with a central peak and some 
broader background structure is present. Also, the 15ms 
period together with the (rotating frame) oscillation at 
frequency A introduces ripples in both the experimen- 
tal and calculated traces with spacing (15ms)" 1 = 67Hz. 
However, there is some asymmetry in the experimental 
trace which naturally is not present in the modeling. This 
asymmetry could be an effect of non-perfect orientation 
or maybe an effect of the frequency scan being a little too 
fast. But this is a minor detail to the general impression 
that the modeling is doing quite well. 

We have tried to vary the duration of the pump pulse, 
the duration of the probe pulse, and the value of r pump . 
In all cases we maintain the agreement between measured 
and calculated spectra as demonstrated for the single ex- 
ample in Fig. El If r pump is strong, the magneto-optical 
resonance signal is essentially a broad feature the struc- 
ture of which in principle can be calculated. But varying 
the total number of atoms we h ave seen exp erimentally 
that a simple relation like J z oc V area • width practically 
holds very well where J z is measured by a DC-Faraday 
rotation signal Q18J1 and the area and width refer to the 
pulsed spectrum. This relation is exact in the cw case 
for p = 1 but is apparently useful in some cases for the 
pulsed spectrum. This is an advantage since we in these 
cases need not perform the more cumbersome modeling 
suggested by JZjJl . 



V. CONCLUSION 

We have studied the magneto-optical resonance for ce- 
sium in a vapor cell as a method for characterization of 
the quantum spin state of the atomic vapor. We have op- 
tically pumped the atoms to various distributions among 
the Zeeman sublevels and exploited the magneto-optical 
resonance signal, in particular the quadratic Zeeman ef- 
fect, to characterize these distributions and the decay 
rate of coherences among the levels. 

The theoretical description relies on well known tech- 
niques combined with ad hoc tailored models for our 
particular needs. Specially useful has been the relation 



Pm,m e m which just happened to be a sufficient model 
in many of our cases. Our experiments agree very well 
with the models in the cw case and demonstrate that 
signals obtained in the pulsed regime are useful and well 
understood. Our approach with the quadratic Zeeman ef- 
fect is particularly well suited for characterization of the 
ground spin state within a hyperfine multiplet. Given 
the off-resonant laser probe with low photon scattering 
rate we obtain the required high resolution. The meth- 
ods could be extended to the microwave region and cover 
transitions between different hyperfine levels similar to 
the approach in [jjj . For excited spin state characteriza- 
tion the fast spontaneous decay rules out high resolution 
methods. In this case scattered photons from the pump- 
ing process may provide useful information [T^ |. 

We have been able to create the coherent spin state to a 
high degree of accuracy (better than 98%) which is a very 
good starting point for studies of quantum effects in our 
spin ensemble. One limitation of our procedure is the fact 
that high orientations and high pumping rates decrease 
the possible resolution as discussed in section llll Dl Thus 
it takes some effort and experimental stability to see the 
difference between a 98% and a 100% polarized sample. 
Also, a pulsed laser setup complicates the conclusions 
but valuable information can still be extracted. The use 
of 894nm pumping light on the £>l-transition turns the 
F = 4, rriF = 4 state into a dark state which has proved 
to be essential to obtain high degrees of orientation. 

The reliable method for characterization of a coherent 
spin state of atoms in a paraffin coated vapor cell de- 
scribed above provides the basis for implementation of 
various entanglement and quantum communication pro- 
tocols utilizing collective atomic spin states. 
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APPENDIX A: THE QUADRATIC ZEEMAN 
EFFECT 

The quadratic Zeeman effect is well understood 0] , we 
will just outline the important results for completeness 
of this paper. An alkaline atom in an external magnetic 
field B is described by the Hamiltonian 

H = hal ■ J - fj-J B — B (Al) 

where J describes the angular momentum of the out- 
ermost electron, I is the nuclear spin, a describes the 
strength of the magnetic dipole interaction between the 
electronic and nuclear spin, and h is Planck's con- 
stant. The magnetic moment of the electron (for an s- 
electron with L = 0) is fij = -1.0011596521869(41)/i B . 
The value for the nuclear moment in cesium is /!/ = 
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2.582025(4)^. Thus, the last term in (IA If) always gives 
a minor correction compared to the second term, but the 
relative strength between the first and second terms de- 
pends on the magnetic field B. 

The exact solution for the energy E to the above 
Hamiltonian is 



Ef,™ — — 



2(21+ 1 



-T Bm± — 



l 



4m 



21 + 1 

(A2) 

where ± is used for F = 1+1/2, m is the magnetic 
quantum number (quantized along the direction of the 
magnetic field), B = |B|, and the hyperfine splitting v^is 
relates to a by hvhfs = 4f (2/ +1). The parameter x 
describes the relative strength between the Zeeman effect 
and the hyperfine splitting: 



{-m/J + fii/I)B 
hu his 



(A3) 



For weak fields m describes the projection of the total 
angular momentum F = I+J. The energy levels (SH can 
be seen in the inset of Fig.^ as a function of the strength 
of the magnetic field. We see that for small field strengths 
or very strong fields, the energy depends linearly on B. In 
the intermediate region the situation is quite non-linear. 
Our experiment is performed in the weak field regime 
with x w 3 • 1CP 4 . Here a linear approximation is very 
good, but it is important to calculate also the second 
order contribution. 

We study the magnetic sublevels by inducing magnetic 
transitions with Am = ±1. Thus, it will be interesting to 
calculate the separation of adjacent sublevels. We start 
out by expanding ilA2j) to first order in the magnetic field 
strength B (leaving out the constant shift independent of 
B). With the standard convention 



E Fm = gFVBBm 



(A4) 



we get for cesium with nuclear spin 1 = 7/2 



9f = 



0.250390 
-0.251194 



-jij/J + hi /I 

21+1 
for F = 4 
for F = 3 



(A5) 



These two numbers differ in magnitude by approximately 
0.3%. Thus, in the lower linear regime of the inset of 
Fig. we have a slightly higher separation between levels 
for the case of F = 3 than F = 4. 

To calculate the quadratic Zeeman shift, it will suffice 
to do the approximation [ii = 0. In this case we may 
write the first order expansion of as hv-L = E m+ i — 

E m = ~2i+i ' an d we then easily derive to second 
order 



E, 



m+l 



E,„ 



VL 1 



(2m + 1) 



(A6) 



This equation describes the transition frequency between 
the m'th and the (m + l)'th level. The separation vqz 
caused by the quadratic Zeeman effect between two 
in e.g. Fig. will thus be 



^QZ 



Vhfs 



(A7) 



All our experiments reported in this paper will have in 
the vicinity of 325kHz corresponding to a magnetic field 
of a little less than 1 Gauss. With the cesium hyper- 
fine splitting being z/hf s = 9.1926GHz we get a quadratic 
Zeeman splitting of 23Hz. 
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